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ABSTRACT 


This  report  presents  a  detailed  description  of  a  computer  program  for  analyzing  crack 
propagation  in  cyclic  loaded  structures.  The  program  calculates  crack  growth,  both  for 
uniform  and  non-uniform  cyclic  loading  and  also  calculates  the  number  of  load  cycles  to 
cause  instability  of  crack  growth.  Instructions  for  use  of  the  program  and  two  illustrative 
analysis  problems  are  presented. 
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SECTION  I 
INTRODUCTION 


The  computer  program  presented  in  this  report  was  developed  to  support  in-house  research 
efforts  by  the  AF  Flight  Dynamics  Laboratory  in  crack  propagation  analysis  of  cyclic  loaded 
structures.  The  program  was  specifically  written  to  take  into  account  complex  crack  geometries 
and  cyclic  loads  of  nonuniform  character.  The  generality  of  the  program  and  the  simplicity 
of  input  data  were  achieved  by  writing  and  running  the  program  under  the  digital  simulation 
program  MIMIC.  The  program  in  its  present  form  was  found  very  useful  in  solving  numerous 
problems  occurring  in  aircraft  type  structures,  such  as  fatigue  crack  propagation  and  crack 
growth  emanating  from  projectile  impact  damage. 
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SECTION  II 

DESCRIPTION  OF  COMPUTER  PROGRAM 


This  program  calculates  the  crack  propagation  behavior  of  cyclic  loaded  structures  by 
means  of  the  theory  described  in  Reference  1.  The  theory  states  that  the  crack  growth  rate 
is  governed  by  the  following  first  order  linear  differential  equation: 


dA  C(DK)SN 

dT  '  (l-R)KC-DK 


For  problems  of  cyclic  load  of  uniform  character,  e.g.  the  load  range,  DP,  and  the  load 
ratio,  R,  constant,  equation  (1)  has  the  form 


dA 

dT 


f  (A) 


(2) 


Given  an  initial  crack  size  at  an  initial  value  of  T,  such  as  T  =  0,  the  computer  program 
calculates  the  crack  length  after  a  given  number  of  cycles.  The  program  also  calculates  the 
number  of  cycles  required  for  crack  growth  instability,  or  the  point  when  the  denominator  in 
equation  (1)  first  becomes  negative. 

The  program  was  written  in  MIMIC,  a  Digital  Simulation  Program  developed  at  Wright- 
Patterson  AFB,  and  was  run  on  an  IBM  7094/7044  DCS  Computer  with  a  Fortran  IV  IBSYS 
monitor.  Thus,  this  program  may  be  run  at  any  installation  which  has  the  capability  of  pro¬ 
cessing  MIMIC  programs.  The  MIMIC  program  is  available  on  request  (see  reference  2)  and 
has  been  written  for  a  number  of  computers. 


Input  to  the  crack  propagation  program  includes  the  following  parameters  for  all  problems  to 
be  analyzed: 


AI  Initial  value  of  A 

C  Material  constant 

SN  Power  of  DK 

KC  Fracture  toughness  parameter 

DTI  Printout  interval  for  T 

TI  Initial  cycle  number 


Other  parameters,  such  as  the  plate  width,  B,  are  usually  required,  but  they  will  depend  on 
the  particular  problem  to  be  solved. 

In  addition,  the  function  representing  the  stress  intensity  factor  range,  DK,  must  be  defined. 
The  function  can  be  expressed  analytically,  or  it  can  be  listed  as  point  values  in  tabular  form. 
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For  example,  many  solutions  for  DK  can  be  expressed  as  follows: 


DK  =  DP  Jtt  A  FX 


(3) 


where  DP  is  the  loading  parameter  and  FX  is  a  correction  factor.  For  the  solution  of  equation 
(1)  ,  DP  and  R  must  always  be  given  in  tabular  listings  as  functions  of  the  variable  T.  The 
correction  factor  FX  can  either  be  given  as  an  analytic  function  of  the  crack  length  A,  or 
listed  as  point  data  in  tabular  form.  If  tabular  form  is  needed,  the  MIMIC  Program  has  the 
capability  of  generating  functions  of  either  one  or  two  variables,  such  as 


FX  =  f ( a ) 
or 

FX  =  f(a,£) 


where 


a  =  a(A),  (3  =  /3(A) 


A  general  description  of  the  MIMIC  Program  with  operating  instructions  is  given  in  Re¬ 
ference  2.  The  card  sequence  for  crack  propagation  analysis  is  given  in  Section  IIL  General 
input  instructions  are  given  in  Section  IV. 
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SECTION  III 
CARD  SEQUENCE 

The  input  to  the  computer  program  consists  of  three  sets  of  cards.  They  are  the  (1)  MIMIC 
system  control  cards,  (2)  function  cards  and  (3)  input  data  cards.  These  are  described  in 
sequence  as  follows: 

1.  MIMIC  System  Control  Cards  -  The  set  of  cards  from  $SETUP  thru  $DELETE. 

2.  Function  Cards  -  These  are  the  program  cards  and  include  the  following: 

a.  Program  Name  Comment  Card 

b.  Constant  Name  Cards 

c.  Tabulated  Function  Definition  Cards* 

d.  Parameter  Name  Card  -  Defines  the  step  functions  DP  and  R 

e.  DK  Definition  Cards 

f.  Differential  Equation  Cards 

g.  Program  Stop  Cards 

h.  Header  and  Output  Cards 

i.  END  Card 

3.  Input  Data  Cards  -  This  set  of  cards  immediately  follows  the  program  END  Card. 

a.  Constant  Cards 

b.  Tabulated  Data  Cards  for  Function  Generators* 

c.  Tabulated  Data  Cards  for  Loading  Parameters 

A  complete  listing  of  the  program  input  for  two  different  problems  is  presented  in  the  Appendix. 


*  Use  only  if  required. 
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SECTION  IV 

GENERAL  INPUT  INSTRUCTIONS 


1.  MIMIC  SYSTEM  CONTROL  CARDS 

The  MIMIC  program  here  at  WPAFB  is  maintained  on  a  tape  (called  MIMIC).  The  control 
cards  call  for  the  mounting  of  this  tape  ($SETUP)  and  for  the  calling  forth  from  this  tape  of 
the  program  subroutines  ($IBLDR).  The  subroutine  MM02  was  modified  for  this  program  to 
change  the  standard  output  format.  This  change  requires  insertion  of  the  modified  Fortran 
source  or  binary  deck  surrounded  by  $IEDIT  control  cards  in  place  of  the  $IBLDR  MM02  card 
in  the  normal  MIMIC  system  control  cards. 

The  control  cards  end  with  a  $DATA  card  to  signal  the  beginning  of  the  function  cards  and 
a  $DELETE  card  to  delete  the  MIMIC  compiled  listing.  A  complete  listing  of  the  cards  is 
shown  in  TABLE  I. 

2.  FUNCTION  CARDS 

a.  Comment  Cards:  Any  symbol  placed  in  colume  1  results  in  the  entire  card  being 
treated  as  a  comment  card.  The  use  of  these  cards  is  optional  in  the  program. 

b.  Constant  Name  Cards:  The  name  of  a  constant  is  defined  by  entering  it  on  a  CON  card, 
and  its  numerical  value  is  given  in  the  data  section  of  the  program.  As  many  as  six  constants 
may  be  named  on  one  CON  card  and  as  many  CON  cards  as  necessary  may  be  used  in  the 


program 

i.  The  format  for  the  Constant  Name  Cards  is  as  follows: 

10 

19 

73 

CON(Al,  C,  SN,  KC,  *,  *) 

CON(*,  *,  *,  *,  *,  *,) 

CON(TI,  DTI,  STOP1) 

The  asterisks  designate  additional  constants  which  may  be  required  for  particular  problems, 
such  as  the  plate  width  B,  or  an  angle  a  .  All  specific  constants  shown  are  required  for  every 
problem. 


c.  Tabulated  Function  Definition  Cards:  Functions  generated  from  tabular  data  are 
specified  by  a  CFN  (constant  function)  card.  The  name  of  the  function  (e.g. ,  F)  and  the  number 
of  pairs  or  triples  of  points,  n,  are  entered  on  the  CFN  cards  as  follows: 


10 

19 

73 

F 

C  FN(n.) 

The  numerical  data  for  the  function  is  given  in  the  data  section  of  the  program.  The  function 
is  used  in  the  program  by  specifying  the  array  name  and  the  independent  variable  on  a  FUN 
card  (See  examples  in  Tables  IV). 


d.  Parameter  Name  Card:  This  card  defines  the  load  parameters  for  the  program.  The 
parameters  are  named  on  a  PAR  card  and  their  numerical  values  entered  on  data  cards  using 
the  same  format  as  for  constants  (See  examples  in  Tables  II  and  IV).  The  basic  parameters 
required  for  every  problem  are  TF,  P,  and  R.  Additional  parameters  can  be  specified  for 
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solving  problems  of  combined  loading,  or  for  problems  where  constants  such  as  C,  SN,  and 
KC  change  due  to  environmental  conditions. 

e.  DK  Definition  Cards:  These  cards  define  the  expression  for  DK  for  the  particular 
problem  to  be  solved.  One  card  should  be  used  to  define  the  expression  for  DK  and  additional 
cards  should  be  used  to  define  the  constants  and  functions  in  the  expression.  The  format  for 
these  cards  is  shown  in  the  following  example  from  Table  IV: 


10 

19 

Result 

Expression 

PI 

3.14159 

K 

(A/B)  *(A/B) 

FX 

FUN(F,  A/B) 

DK 

DP*SQR(PI*A*FX*FX) 

f.  Differential  Equation  Cards: 


g.  Program  Stop  Cards: 

h.  Header  and  Output  Cards: 

i.  END  Card: 

3.  INPUT  DATA  CARDS 

Input  format  for  these  cards  is  given  in  Reference  2.  Examples  for  two  different  problems 
are  shown  in  Tables  II  and  IV.  As  the  examples  indicate,  the  input  requirements  for  these 
cards  are  that  numerical  values  must  be  entered  in  the  same  order  as  the  constants  or 
parameters  on  the  CON,  CF,  or  PAR  cards.  The  first  number  must  be  entered  in  columns 
1-12,  the  second  in  columns  13-24,  the  third  in  columns  25-36,  the  fourth  in  columns  37-48,  the 
fifth  in  columns  49-60,  and  the  sixth  in  columns  61-72.  Finally,  the  numbers  must  be  written  in 
floating  point  form,  that  is,  a  number  of  digits  with  a  decimal  point  somewhere  in  the  number. 


These  cards  do  not 
change.  See  examples 
in  Tables  II  and  IV. 
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APPENDIX 

1.  Sample  Problem  of  Griffith  Crack  in  a  Finite  Width  Plate 

a.  Problem  Description:  See  Figure  1. 

b.  Load  Spectrum:  See  Figure  3. 

c.  Listing  of  Data  Input:  See  Table  II. 

d.  Output  of  Results:  See  Table  III. 

2.  Sample  Problem  of  Crack  Emanating  from  a  Circular  Hole 

a.  Problem  Description:  See  Figure  2. 

b.  Load  Spectrum:  See  Figure  3. 

c.  Listing  of  Data  Input:  See  Table  IV. 

d.  Output  of  Results:  See  Table  V. 


9 


AFFDL-TR-67-5 


DK  =  DP^/tt  AyS  2  (from  Reference  I) 

where 

/32  =  I  +  I.  I8968K2  +  I  30I62K4  +  I.36502K6 

+  I  37394K®  +  I  47638K10 

and 

K  =  A/B 


Figure  1.  Griffith  Crack  in  a  Finite  Width  Plate 
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DK  =  OPjTT  FX  (From  Reference  3) 


where  point  volues  of  FX  are  as  follows. 


A/B 

FX 

0  0 

3  39 

0  10 

2  73 

0.20 

2.30 

0.30 

2  .04 

0.40 

1  86 

0  50 

1  73 

0.60 

1.64 

0.80 

1.47 

1.0 

1  37 

1.5 

1  18 

2.0 

1  06 

3  0 

0.94 

5.0 

0.81 

10.0 

0  75 

00 

0  707 

Figure  2.  Crack  Emanating  from  a  Circular  Hole  in  a  Sheet 
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Figure  3.  Description  of  Cyclic  Loads  for  a  Sample  Problem 
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TABLE  I 

LISTING  OF  INPUT  CARDS 


SSJETUP_J_B4  MIMIC  »NORI  NG 

SASSIGN  SYSLB4 


SI B JOB 

MIMIC  F IOCS 

SI  ED  I  T 

SYSLB4, SCHF1 

SIBLDR 

MMC  I 

SI  EDIT 

SI BFTC 

MMO  2  XR7 

SUBROUTINE  MIMEX 

MM0200 1 0 

C 

MM020020 

EXECUTION  PROGRAM 

MMC2r'03° 

C 

MMO20^4n 

DIMENSION  P( 95 ) ,R( 250C > , S ( 2500 ) t BCD( IC , 900 ) , FF ( 9 1^0 ) 

MM02r'°5n 

COMMON  P,R,S,FF,IOUT,IPAR,  INOUT,  IHDR, I F I N , I  END , NP AR , 

MM02^06n 

c 

:  DUMMY(3)*IPC»IS»  DUMMY 1 

MM02DO70 

EQUIVALENCE  ( BCD,  FF ( 10  1 )  ) 

MMO200R0 

NRUN-0 

MM020090 

I  NOUT  =99 

MMH20 ion 

CALL  MIMIO(A,B,C»D,E»G) 

MM020I 1 o 

c 

MM020I 20 

c**** 

PROGRAM  EXECUTION  CONTROL  SWITCHES 

MMO  2n I 30 

c 

MM020140 

c 

I END=0  NOT  END  OF  RUN 

MMO  20 1 50 

c 

=L  END  OF  RUN 

MMO  20160 

c 

I P A R  =  0  DU  NOT  READ  PARAMETERS 

MMO 20 1 70 

c 

=  1  READ  PARAMETERS 

MMO20I 80 

c 

I OUT  =0  DO  NOT  WRITE  OUTPUT 

MM0201 90 

c 

=i  WRITE  OUTPUT 

MMO20200 

c 

I F I N=C  DU  NOT  TEST  FINISH  STATEMENTS 

MM020210 

c 

=1  TEST  FINISH  STATEMENTS 

MM020220 

c 

I HDR=0  DO  NOT  WRITE  HEADERS 

MM02023O 

c 

=L  WRITE  HEADERS 

MM02024O 

c 

NP AR=C  NO  PARAMETERS  (ONF  RUN) 

MMO20250 

c 

=1  PARAMETERS  (ONE  OR  MORE  RUNS) 

MM02026O 

c 

MMO?027^ 

SET  FOR( READING  INPUT  DATA 

MM0202BO 

c 

MM02  0290 

IOC 

NRUN=NRUN+ I 

MMO203OO 

I  P AR  =  I 

MM0203IO 

PRINT  I  * NRUN 

MMO20320 

I END=0 

MMO20330 

I  OUT  =0 

MMO20340 

I F I N=0 

MM020350 

I H  D  R  =  0 

MM02O360 

R  (  I  )  =  0  . 

MM02O37O 

c 

WR I T  E ( 6  t  2 ) 

MMO2039A 

CALL  F 

MMQ20390 

i par=d 

MMO204OA 

CALL  F 

MM'>20410 

SET  FUR  HEADING,  OUTPUT,  TESTING  END  OF  RUN 

MM02042O 
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TABLE  I  (Continued) 
LISTING  OF  INPUT  CARDS 

IHOR= 1 
200  I F I N= 1 
I OUT  = 1 
CALL  F 
I HDR=0 
I F I N=0 
I 0UT=0 

C****  JEST  FOR  END  OF  RUN 

I F ( IEND.NE.0  )  GO  TO  300 
C 

C****  INTEGRATE 
C 

CALL  MIMIN 
GO  TO  200 
C 

C***«  TEST  FOR  FURTHER  PUNS 
C 

300  IFdPC.EQ.il  CALL  MIMLT 
IF(NPAR.NE.C)  GO  TO  ICO 
RETURN 

1  FORMAT ( 1 1H  BEGIN  RUN  ,14) 

2  FORMAT { 1H1 ) 

END 

SIEDIT 
$  I BLDR 
SIBLDR 
$1 BLDR 
$  I BLDR 
$  I BLCR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SIBLDR 
SI  ED  I  T 
$DATA 

I 

Data  Input  (See  Tables  II  and  IV) 

{ 

SEOF 


SYSLB4, SCHF 1 

MM03 
MM04 
MM05 
MMC6 
MMG7 
MM08 
MM^9 
MM1C 
MM1 1 
MM  1 2 
MM  13 
MM14 
MM  1 5 
MM16 
MM  17 
MM  1 B 
MM19 
MM95 
MM96 
MM97 
MM9B 
MM99 


MM02O43O 
MM02n44O 
MM02045r 
MM92046O 
MM020470 
MMO204B0 
MM02O49'-' 
MM020  500 
MM02C  51 0 
MM02052O 
MM020530 
MM020540 
MM02C550 
MM020560 
MM020570 
MMO  20  5  80 
MM°2r'590 
MM020600 
MM02061Q 
MM020620 
MM020630 
MM020649 
MM020650 
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TABLE  II 

SAMPLE  PROBLEM  OF  GRIFFITH  CRACK  IN  A  FINITE  WIDTH  PLATE,  LISTING  OF 
DATA  INPUT 


S  DAT  A 

C  GRIFFITH  CRACK  IN  FINITE  WIDTH  PLATE 
SDELETE 

CON ( A  I »  C  »  SN  *  KC * B ) 

CON ( T I  »  DT 1 »  STOP  1 ) 

C  hTOPl  IS  A  CONTROL  CONSTANT  WHICH  WILL  STOP  THE  PROGRAM  AFTER 
C  DENOM  HAS  BECOME  NEGATIVE. 

PARI TF,DP»R ) 

PI  3.14159 

DTMIN  1.0 

K  ( A/B ) * ( A/B ) 

BETA1  1 .0+1. 18968*K+1.30162*K*K 

BETA2  K*K*K*( 1.36502+1. 37394*K+ 1.4763 8* K*K) 

B2  BETA1  +  BETA2 

DK  DP*SQR(PI*ABS(A)*B2) 

NUMER  EXP  I SN»DK*EXP I  1 . /SN »C )) 

DENOM  ( 1 . -R ) *KC-DK 


I  DA 
A 

DT 

STOP  ST  OP  1 

STOP 

STOP  A I 

STOP  T I 

LCV 

LCV 

LCV 


NUMFR/DENOM 
INTI  I  DA  » A  I ) 

FSW I TI+T  +  DT1-TF ,DT1 ,DT 1 , TF-T  I -T ) 

FSW I DENOM  »-l • 0  *  — 1 • 0  »STOP 1 ) 

IOR(FIN(TI+T,TF),FIN(O.O.STOP1)  » F I N ( 0 . 0  * DENOM )  ) 
A 

TF 

FSW I  STOP  1  ♦  FALSE .FALSE,  TRUE) 

HDR I T, A, R, DP. DENOM) 

OUT(TI+T»A»R»DP»  DENOM ) 

END 


1 .25 
0.0 
1.0 

1 ooo.n 
1 001 .0 
POOO.O 
*EOF 


5.E-13 

1000. 

6000.0 

2000.0 

4000.0 

4000.0 


3.0 
1.0 
0 .0 
0.666 
0.500 
0.333 


68000.  10. 
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TAELE  IV 

SAMPLE  PROBLEM  OF  CRACK  EMANATING  FROM  A  CIRCULAR  HOLE,  LISTING  OF  DATA  INFUT 


AT  A 

C  CRACK  EMANATING  FROM  A  CIRCULAR  HOLE. 

^DELETE 

CON ( A I *C»SN*KC*B* ) 

CON ( T I  *  DT 1 »  STOP  1 ) 

C  STOP1  IS  A  CONTROL  CONSTANT  WHICH  WILL  STOP  THE  PROGRAM  AFTER 

C  DENOM  HAS  BECOME  NEGATIVE. 

F  CFN ( 15.  ) 

C  F  IS  GIVEN  IN  TABULAR  FORM. 

PAR ( TF »  DP  » R ) 

PI  3.14159 

DTMIN  1.0 

K  (A/B)*(A/B) 

FX  FUN ( F »  A/ B ) 

DK  DP*SQR(PI*A*FX*FX) 

NUMER  EXP(SN»DK*EXP(1./SN»C)) 

DENOM  (l.-R)*KC-DK 


STOP 

STOP 

STOP 

LCV 

LCV 


IDA  NUMER/DENOM 

A  I  NT (  I  DA  » A  I  ) 

DT  FSW( TI+T+DT1-TF ,DT1»DT1»TF-TI-T) 

STOP  1  FSW ( DENOM  » - 1 • 0  » - 1 • 0  »  STOP  1 ) 

STOP  I  OR (FIN(TI  +  T»TF) » F I N < 0 . 0 » ST0P1 )  ,  F I N ( 0 . 0  * DENOM )  ) 

A I  A 

T I  TF 

LCV  FSW( ST0P1 »FALSE»FALSE .TRUE ) 

HDR ( T, A, R, DP, DENOM) 

OUT( TI  +  T»A,R *  DP, DENOM) 

END 


1.0 

5.0E-13 

3.0 

0.0 

1000. 

1.0 

0.0 

3.39 

0.10 

2.73 

0.20 

2.30 

0.30 

2.04 

0.40 

1.86 

0.50 

1.73 

0.60 

1.64 

0.80 

1.47 

1.0 

1.37 

1.5 

1.18 

2.0 

1.06 

3.0 

0.94 

5.0 

0.81 

10.0 

0.75 

100.0 

0.707 

1.0 

6000.0 

0.0 

1000.0 

2000.0 

0.666 

1001.0 

4000.0 

0.500 

8000.0 

4000.0 

0.333 

SEOF 

66000.0  0.25 
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STOP  T I  TF 

LCV  FSW(ST0P1, FALSE, FALSE, TRUE) 
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